Disorder effects in the AHE induced by Berry curvature 
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We describe the charge transport in ferromagnets with spin orbit coupled Bloch bands by combin- 
ing the wave-packet evolution equations with the classical Boltzmann equation. This approach can 
be justified in the limit of a smooth disorder potential. Besides the skew scattering contribution, 
we demonstrate how other effects of disorder appear which are closely linked to the Berry curvature 
of the Bloch states associated with the wavepacket. We show that, although being of the same 
order of magnitude as the clean limit contribution, generally disorder corrections depend differently 
on various parameters and can lead to the sign reversal of the Hall current as the function of the 
chemical potential in systems with a non-constant Berry curvature in momentum space. Earlier 
conclusions on the effects of disorder on the anomalous Hall effect depended stricly on the lack of 
momentum dependence of the Berry curvature in the models studied and generalizations of their 
findings to other systems with more complicated band structures were unjustified. 

PACS numbers: 73.20.Fz, 72.10.Fk,72.15.Rn 



I. INTRODUCTION 

The theory of the anomalous Hall Effect (AHE) has a 
long history The appearance of the Hall current requires 
breaking of some basic symmetries. It was proposed by 
Karplus and Luttinger in early 50si that the anomalous 
Hall effect in ferromagnets results from the interplay of 
the exchange field, which breaks the time-reversal sym- 
metry, and the spin-orbit coupling, that violates the chi- 
ral symmetry. Interestingly, at the same time a similar 
effect was predicted and explained in geometrical optics 2 , 
however, its relation to the AHE was revealed only 
recently^ after the modern interpretation of both effects 
in terms of the Berry phase, which affects the motion of 
wave packets, had been constructed iMii Luttinger— built 
a detailed theory of the AHE based on the high order 
quantum Boltzmann equation calculations— In that work 
Luttinger identified various contributions, known today 
as Berry phase contribution, the skew scattering, the side 
jump on the impurity potential and a contribution that 
involves interference from many scatterers. 

Since then, a number of theoretical works appeared 
that extended the theory. However untill recent time 
most of them had been devoted to what is today called 
the extrinsic AHE. In the extrinsic AHE a simple Bloch 
band structure of the system is assumed and the spin or- 
bit interaction is localized on the impurity potential via 
terms of the form Aso^zk x VVdis(x) e.t.c. It was rec- 
ognized by Smit, Berger and otherc^VWW 5 ! 16 that 
in this case the main contributions to the Hall current 
will be those from impurities via the side jump and the 
skew scattering mechanisms. In contrast to the extrin- 
sic AHE, the intrinsic one assumes that the spin orbit 
coupling is already present in the band structure of the 
system and generally cannot be considered as weak in 
comparison even with the Fermi energy. Recently the 
interest toward the intrinsic AHE has grown up consid- 
erably due to new applications in the diluted magnetic 
semiconductors (DMS) and due to the interesting modern 



interpretation of the Hall current in terms of the Berry 
phase— iii The Berry phase contribution to the Hall cur- 
rent was shown to be in good quantitative agreement 
with experiment in many different materials with strong 
spin-orbit couping in their band structure, giving weight 
to the theory involving this contribution alone— Li&iS 

The Luttinger's theory of the intrinsic AHE, however, 
predicts that other contributions should arise due to the 
scattering from a disordered potential even if the disor- 
der potential itself is spin independent. The asymmetry 
in scatterings becomes transparent in the basis related 
to the Bloch states. As an example, Luttinger consid- 
ered a rather simplified model which demonstrated that 
such corrections must reverse the sign of the Hall current 
in the dc- limit, in comparison to the clean limit, though 
in the high frequency AC-case the clean limit contribu- 
tion should dominate. The same results have been shown 
by several other works utilizing different approaches and 
focusing on this simplified model— 

The results from these simplified models seem to be in 
a contradiction with recent work which did not find such 
a change of sign in numerical simulations or in a compar- 
ison with the experimental results. We show here that a 
possible resolution of this discrepancy is the simplicity of 
the early models from which many generalizations were 
stated without justification and whose results depended 
drastically on the simple momentum depencence of the 
Berry curvature of the Bloch states. 

The work by Luttinger and works of other authors, 
related to the disorder contribution to the AHE, are 
rather involved. Generally various contributions sepa- 
rately turned out to be not gauge invariant and only the 
final result was physically meaningful. Because of these 
shortcomings, we reformulate the basic arguments of the 
previous authors in terms of wave packet dyanamic equa- 
tions which are fully gauge invariant and consistent with 
prior results*! Being gauge invariant, the wave packet 
equations allow one to identify the physical meaning of 
various contributions. We should note, however, that 
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the wave packet equations are valid only in the limit of 
smooth potentials. This restricts the applicability of our 
conclusions and generalizations of our results should not 
be made to regimes where the semi-classical treatment is 
not justified without careful checks; therefore we did not 
make the goal to construct the final theory of the disor- 
der in the AHE but rather to construct a simple formal- 
ism that demonstrates basic features of the problem and 
highlights a basic key ingredient missed by prior theories 
studying the effects of disorder, namely the importance 
of the momentum dependence of the Berry curvature. 

We have found that the change of the sign of the Hall 
current due to the disorder, found in the early works, is 
not universal and arises in the models where bands have 
a constant Berry curvature in the momentum space. In 
diluted magnetic semiconductors the Berry curvature is 
strongly momentum dependent. This separates paramet- 
rically the clean contribution from the others and makes 
the sign reversal not universal. Rather generally we find 
the same sign of the total Hall current for realistic bands 
except in extreme situations where one of the bands be- 
comes depleted. 

We organize the rest of the paper as follows. In Sec. II 
we review the wave packet dynamics theory of Sundaram 
and Niu 6 . In Sec. Ill we analyze the model for constant 
Berry curvature obtaining in a physically clear way the 
reversal of the sign in the presence of smooth disorder. 
In Sec. IV we apply the theory to the case of the Rashba 
model and show the non-trivial dependence of the con- 
tributions from disorder scattering and the clean Berry 
phase contribution as a function of the Fermi energy and 
in Sec. V we present our conclusions. 

II. THE WAVE PACKET EQUATIONS 

The motion of wave packets formed by Bloch states is 
governed by the following equations. 6 

|k = eE-W(r) (1) 

d de(k) dk „ 

— r = — — x F (2) 

dt dk dt y ' 

where e(k) is the energy dispersion in the band, E is 
the external electric field acting on a wave packet having 
the electric charge e, V(r) is the local potential in the 
sample, for example the potential of impurities and F 
is the Berry curvature of the Bloch band. For the two- 
dimentional motion only the out of plane component of 
F is nonzero 



The second term in (|2J is responsible for the so called 
anomalous velocity. For the reader not familiar with the 
notion of the anomalous velocity we provide the appendix 
with a simple example from classical physics that gives 
an intuitive explanation of the physical meaning of the 
anomalous velocity in Rashba coupled 2DEG. For the 
rigorous theory we refer to the original papers^. 

The anomalous velocity is orthogonal to the direction 
of the electric fieldi^ (which we chose to be along the 
x-axes) . 



The semiclassical equations (QJ and © map the quan- 
tum mechanical problem to a classical one where par- 
ticles have the electric charge e and move according to 
equations of the wave packet dynamics. We assume that 
at equilibrium the distribution of such classical particles 
is the same as the Dirac distribution of electrons in the 
sample. This mapping to a classical system considerably 
simplifies the treatment of the problem both analytically 
and numerically. Analytically, one can apply the classical 
Boltzmann equation approach to calculate the transport 
coefficients, numerically the molecular dynamics simula- 
tion of the motion of classical particles is simple and may 
not be restricted to a small system size as, for example, 
in the case of a numerical diagonalization of a quantum 
mechanical Hamiltonian. 

Unfortunately equations and @ are valid in the 
adiabatic limit only. They cannot be applied to the case 
of scatterings on a short range delta- function like po- 
tential. Only scatterings on impurities, whose potential 
varies appreciably only on distances much larger than the 
size of the wave packet can be calculated this way 6 ". In 
spite of this restriction of its applicability the limit of a 
smooth potential is a very interesting one to investigate 
the influence of disorder on the anomalous Hall effect. In 
realistic applications a system with long range impurities 
is realized, for example, in the high mobility 2D electron 
gas 2 * and the out of plane Zeeman field can be induced 
there by polarizing nuclear spins or by introducing addi- 
tional magnetic impurities. 

Recently another related effect, namely the intrinsic 
spin Hall effect, was introduced^. A number of theoreti- 
cal papers have explored the importance of the disorder. 
The debates on this topic are ongoing (see e.gi 26 i 27 i 28 ). 
Understanding the AHE may shed light on the disorder 
role in the spin Hall effect. 



III. CONSTANT BERRY CURVATURE 



A. Clean limit 

(3) 

In many recent applications the Berry curvature F 
here \u s ) is the Bloch state in the absence of the electric strongly depends on the momentum of the wave packet 
field and impurities and s is the index of the band. ^.^1^30 However, it is instructive to consider first the one 
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Bloch band 2D system with a constant Berry curvature 
F z . 

For a constant F z the anomalous velocity (@J leads to 
the following contribution to the Hall current from the 
unperturbed part of the distribution function in a single 
band at zero temperature: 

jiclean) = & J foft^W = 

f kF f 2 * d4> e 2 E x F z kl , . 

= e / kdk / 7^Y2 eE * F * = 5 

k Jo (2^) 2 4^ 

where /o(k) is the equilibrium distribution function and 
kp is the Fermi momentum. The upper index (clean) 
means that the quantity is calculated in the absence of 
disorder and the distribution function of the wave packets 
coincides with the one in equilibrium at E x = 0. 

The case of a constant F z is realized in the conduct- 
ing bands whose states are weakly hybridized with the 
states in the valent hole bands^i where the spin orbit 
coupling is allowed. If the kinetic energies of electrons 
in the conducting band are much smaller than the gap 
between conducting and valent zones, the Berry curva- 
ture due to such an induced spin orbit coupling can be 
considered as a constant for all conducting electrons. In 
modern applications it can be realized in some limiting 
situations, like in the case of R2DEG at a strong Zceman 
field (see following sections). In this section we show 
that for smooth impurity potentials the wave packet ap- 
proach easily reproduces the sign reversal of the Hall cur- 
rent when F z — const. In comparison to previous works, 
however, our approach keeps the derivation gauge invari- 
ant and hence is physically clear. 

B. Effects of disorder 

According to JSJ if the distribution function is the same 
as in the equilibrium at zero external field a Hall cur- 
rent appears in the external field due to the anomalous 
velocity. However, in the steady state the distribution 
function is no longer /o(k). In the electric field and on 
time scales much larger than the scattering time particles 
diffuse with a constant velocity rather than accelerate as 
in the absolutely clean case. The anomalous velocity and 
hence the clean contribution to the Hall current are pro- 
portional to the acceleration k x . Since in the steady state 
the average acceleration is zero up to the first order in 
external electric field E x one can expect that disorder 
should strongly influence the Hall current. 

A natural way to study the effect of disorder on the 
transport in our case is the classical Boltzmann equa- 
tion. There is however a complication when both nonzero 
Berry curvature and finite sizes of impurities must be 
considered. At a scattering on an impurity potential not 
only the momentum but also the coordinate of a particle 
changes. Usually such a coordinate shift at the scatter- 
ing is discarded since after averaging over many scatter- 
ings such random shifts cancel each other and only the 



changes of the momentum matters. However when the 
Berry curvature is nonzero the additional (anomalous) 
shift does not disappear after the averaging. To see this, 
suppose that the term with the Berry curvature F z in 
Eq. J2J) is small in comparison with the first one and 
calculate the corresponding correction to the shift of the 
particle during the scattering on an impurity. Integrat- 
ing |J2J over the time interval at which a particle feels the 
impurity potential during a single scattering and treat- 
ing the second term in (J2J as a small perturbation one 
can find that, after the scattering, a particle makes an 
additional shift 

Ar(k, k')=zx (k - k)F z + J J ^dt (6) 

here k' and k are momentums respectively after and be- 
fore the scattering, t\ and ti are times of entering and 
leaving the impurity in a semiclassical picture. The sec- 
ond term in the rhs. of © is just the shift due to the 
normal velocity. To first order in F z and W it is not 
affected by the Berry phase and hence is averaged to zero 
after many scatterings; therefore we will disregard it in 
our future discussion. 

The particle's displacement due to the first term in © 
is due to the anomalous velocity. This shift does not 
depend explicitly on the details of the impurity potential 
and does not have the chiral symmetry. There are two 
main rather distinct effects due to the appearance of this 
anomalous shift. The first effect, the so called side-jump, 
is that the y-component of this shift does not cancel after 
the averaging over many scatterings and thus contributes 
to the drift velocity perpendicular to the electric field, 
i.e. to the Hall current. We will focus on this effect in 
the following subsection. The second effect is that when 
a scattering takes place in the presence of an external 
electric field there is a change in the potential energy 
upon a scattering given by 

AU = -eE x Ax (7) 

where Ax is the shift along the external electric 
field ilJJ^LiL Both of these effects, as shown below and 
in the next subsection, give the same contribution as Eq. 
(JSJ but with an opposite sign in the particular case of a 
momentum-independent Berry curvature. Because both 
contributions are linear in E x we are able to consider 
them separetly when considering the linear response of 
the system. 

Let us focus first on the second effect. Here only the 
anomalous part of this shift is important for Hall current 
calculations 

Ax = -F z (k' y - k y ) (8) 

and the effect of the normal part of the shift only renor- 
malizes the diagonal current response to E x . Since the 
total energy remains the same after an elastic scattering, 
the kinetic energy must change by the amount 

Ae(k, k') = e(k') - e(k) = eE x Ax (9) 
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This effect leads to an instability of the initial equilib- 
rium distribution function /o(k) = /o( e (k)) due to scat- 
terings. The Boltzmann equation is given by 

% = - E ^ ^ k )) - / ^ k '))] ( 10 ) 

where w(k, k') is the scattering rate. Kinetic energies 
before and after a scattering do not coincide but the dif- 
ference between them is small. Since we are seeking a 
contribution linear in E x we can substitute / by /o which 
depends only on the kinetic energy e(k) so we can expand 



/o(e(k))-/ (e(k')) = 



#(e(k')-e(k)) 



-^eE x F z (k y k' y 



(11) 



This yields 

at 



^(k,k')-^eEF z (k y -k' y )^0 (12) 



The distribution function will relax until a contribution 
to it compensates this relaxation in the stationary limit. 

To find the new equilibrium distribution one should 
substitute /(k) = /o(e(k)) + 5b (k) into the rhs. of equa- 
tion I|1UI) . To make df/dt = 0m equation l|12l) the cor- 
rection 5_e(k) must be the following 



9e (k) = ^-eE x k y F z 



(13) 



This contribution is not symmetric in the y-direction. 
This means that already the normal (i.e. usual <9e/<9k) 
velocity can contribute to the Hall current in the station- 
ary state: 



jnorrnal 



d 2 k 

e 2 E x F z 

(47T) 



g E {K){v y normal) ) 



hp 



(14) 



which has the opposite sign and is exactly the same in 
the absolute magnitude as the anomalous velocity con- 
tribution in the clean limit given by Eq.JSJ). 



C. The side jump 

The anomalous change of energy after the scattering is 
not the only important effect of the anomalous shift J7J. 
The anomalous shift during the scatterings has generally 
a component perpendicular to the direction of the electric 
field. 



Ay(k,k') = F z (k' x - k x ) 



(15) 



It does not contribute to the change of energy during the 
scattering but it shifts a particle along the y-axes. If such 



shifts do not compensate each other after many scatter- 
ings, they should contribute to the total Hall current. 
This phenomenon is known in the theory of the extrinsic 
Hall effect as the side jump. If scatterings happen with 
the rate w(k, k') in average the particle moving with the 
momentum k also acquires the anomalous drift velocity 
perpendicular to the electric field 



(k)\ = <"(k, k')Ay(k, k') = -F z k x /r\ 



where 



II 

(16) 



l/r||=5>(k,k')(l-cos(k,k')) (17) 



and " (sj) n marks the side jump contribution to the phys- 
ical quantity. 

In the equilibrium the side jump does not lead to the 
Hall current because the anomalous velocity (|16fl changes 
the sign under the transformation k x — > —k x . and the 
distribution function in the equilibrium is invariant under 
this transformation. However, when the electric field is 
applied the non-equilibrium correction to the distribution 
function appears which has no such a symmetry under 
the momentum reflection. 

This correction to the distribution function can be de- 
rived by means of the standard approach to the Boltz- 
mann equation. Up to the first order in the electric field 
the standard Boltzmann equation for one band, ignoring 
the asymmetric contribution considered in the previous 
subsection, is 

-eE x v£° rmal) {-^) = -^(k,k')(/(k)-/(k')) 

(18) 

where /o(k) is the equilibrium distribution function. 
Equation 118fl has a solution /(k) = /o(k) + g(k) where 
to the first order in the electric field 



5(k)= e ^T||(-^)4~ J ) 



(19) 



^tjiorma j . g ^ e norma i velocity along the electric field, 
which in our case is 



^normal) = = ^ 



dk x 



(20) 



here we assumed that the conducting band has a trivial 
dispersion e(k) = bk 2 /2 where b = l/m e . The side-jump 
contribution of impurities to the anomalous velocity leads 
to the Hall current 



e 2 E„.Fr f 27r 



(2tt)2 



e 2 E x F., 



x± z 1.2 



(47T) 



(2n) 2 

cos 2 ((/)) J de k k 2 S(e F - e k ) 
(21) 
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As in the previous impurity contribution, the Hall cur- 
rent due to the side jump has the opposite sign and the 
same magnitude as the one in the clean limit. The side 
jump current appears because the anomalous velocity is 
proportional to the acceleration of the particle but in the 
stationary state the average acceleration should be zero 
(up to the terms of the first order in the electric field), 
namely, while the external electric field accelerates the 
wave packet during its motion between impurities, the 
impurity potential generally decelerates it. During such a 
deceleration the wave packet has the anomalous velocity 
with the opposite sign and hence the side jump contri- 
bution should compensate or at least decrease the pure 
limit result. The described picture is of course oversim- 
plified, which will be clear from the discussion of a mo- 
mentum dependent F z . In reality, the external electric 
field equally accelerates all electrons in the Fermi sea but 
effective deceleration due to impurity scatterings is, in 
average, seen only by the electrons near the Fermi level 
because there the distribution function acquires a cor- 
rection required to compensate the acceleration by the 
external field. Hence the side jump contribution is the 
property of the Fermi surface while the clean limit Berry 
phase contribution appears from all electrons even deep 
in the Fermi sea. 



D. Total current and numerical check 

The total Hall current (not including effects of skew 
scatterings due to asymmetry of the collision term kernel) 
is the sum of the clean limit contribution, of the side 
jump and of the normal velocity contributions. For the 
constant F z we find 

j(total) j(clean) _■_ j(sj) _■_ j(?iormal) _ rj x r z ,2 

J yx — J yx J yx ' J yx — K F 

(22) 

As it was predicted in the former literature which focused 
on models with a constant F z , the total current has the 
same magnitude but with the opposite sign as the one in 
the absolutely clean system. Our derivation however is 
considerably simpler and is straightforward to generalize 
to a case with a more complicated band structure such 
as Rashba coupled 2DEG. 

To confirm the analytical result (1221) numerically we 
simulate the motion of particles according to equations 
of motions and Initially all particles were pre- 
pared uniformly distributed over the momentum space 
having the absolute value of the momentum less than a 
specified Fermi momentum. The action of an impurity 
was simulated by the potential of a wall of a finite thick- 
ness with a very fast linearly growing potential inside the 
wall. The wall is sufficiently long to disregard the effects 
on its edges. This type of impurity does not allow the 
skew scattering mechanism to appear because for every 
scattering on the wall that changes momentum from k to 
k' there is a process when a particle hits the same wall 




FIG. 1: Scattering of a particle on impurity. The impurity 
is assumed to have the shape of a long stripe with linearly 
growing potential inside. 

from another side and scatters from k' into k. Thus both 
processes have equal probability and w(k, k') = aj(k', k). 
Figffl demonstrates a scattering on such an impurity. We 
would like to point to the analogy of the scattering in 
Fig^with Goos-Hanchen's and with Fedorov's shifts in 
geometrical optics^. 

One can consider the effect of the wall on a wave 
packet as the one due to an electric field E± acting on it 
only inside the wall in the direction perpendicular to the 
wall. Although the potential in the wall can be strong in 
comparison with the external electric field, we suppose 
that the motion inside the wall is still governed by the 
equations of the wave packet dynamics. Such a scatter- 
ing problem can be solved exactly. This solution shows 
that for the given type of scatterings the theory of the 
anomalous shift is valid even if we do not assume that 
the anomalous velocity term in wave packet equations is 
small. This is useful because the bigger anomalous ve- 
locity is needed to accelerate the numerical calculations. 

Let components of the momentum of the particle in- 
cident on the wall be fc|| and k± parallel and perpen- 
dicular to the wall respectively. In the presence of an 
additional external electric field the total force has gen- 
erally also a nonzero component parallel to the wall eE\ | 
where E\\ is the projection of the external electric field 
on the direction along the wall. Expressions for the fi- 
nal momentum and coordinates right after the scattering 
are strongly simplified in the limit of large E±, so that 
all terms 0(1/E±) can be dropped. In this limit the 
scattering time is vanishing and in @ one can disregard 
the first term in comparison with the second one since 
k± ~ E± is large. Dropping the term with de/dk± the 
equation @ can be readily integrated over the time of the 
scattering leading to the relations Ar|| = eF z (k± — k' ± ) 
and fcn = fc| | . Solving them together with the energy 
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conservation equation e(k) + eE\\r\\ = e(k') we arrive at 
following expressions for the momentum of the outgoing 
particle 



k '\\ = k \\ 

fej. = -k ± -EnF z 



(23) 



In addition, the scattered particle appears in a point of 
the interface shifted in comparison to the incident point 
by the amount. 



An 



2k,F,+EnF 2 



(24) 



The shift r\ | is exactly the anomalous shift in the impurity 
potential. Note also that, as it follows from (|23p. when 
a scattering happens in the additional external electric 
field with nonzero projection to the direction of this shift 
En, the absolute magnitude of the momentum changes, 
which leads to the instability of the equilibrium distribu- 
tion function, discussed in previous sections. The scat- 
tering time in this limit is small enough to be disregarded 
so numerically such a scattering can be easily simulated 
by choosing the coordinate system related to the wall ori- 
entation and updating the momentums and coordinates 
according to (fey, fcj_, r\\, r±) — > (k',, ,kj_,r\\ + Ar||, r_i_), r± 
is the same for the incident and the out going particle. In 
between scatterings equations of free motion in the ex- 
ternal electric field are also trivially solvable. In our sim- 
ulations we assumed the " noncrossing approximation" , 
namely we suppose that every scattering happens on a 
different impurity. The algorithm consists of two parts. 
First we generate randomly the distance L to the next 
impurity. We suppose that between impurities particles 
move only under the action of the electric field. The in- 
tegration of wave packet equations leads to the following 
trajectories 



k x — k x ~\~ E x t 



y .«y 

Ext 2 

y + (k y + F z E x )t 



(25) 



The time t of the motion to the next impurity was esti- 
mated from the equation ((x' — x) 2 + (y' — y) 2 )) 1 ^ 2 = L. 
Solving this relation with the solution (|2*5|l we find the 
momentum and coordinates of the particle before the new 
scattering. In the second step we randomly generate the 
orientation angle of the impurity, switch into the related 
coordinate system and update particle's coordinates and 
the momentum according to l|23|) and l|24|) . Then we re- 
turn to the initial coordinate system and repeat the cir- 
cle. We repeat this circle sufficiently many times to allow 
the distribution to relax to the stationary one (usually a 
few scatterings is enough), however the total evolution 
time was small enough in order to avoid strong heating 
of the system. In our simulations, every particle makes 
about 10 2 — 10 3 scatterings during the whole time. To 
prevent strong heating the electric field is chosen suffi- 
ciently small E x l/Ep ~ 10~ 4 where I is the typical scat- 
tering length and Ep is the Fermi energy. The total 




FIG. 2: Hall conductivity vs strength of the Berry curvature 
in the model with momentum independent F z . Solid line is 
the theoretical prediction that includes impurity scatterings 
and dotted line is the prediction of the clean limit contri- 
bution. Triangles show results of the numerical simulations. 
Simulations were performed in units e = h — m — 1 and 
k F = 1. 



current can be derived as the sum of total displacements 
of all particles divided by the evolution time. 

In FigfJ we compare the analytical result 11221) with 
our numerical simulations for the fixed Fermi energy but 
different strengths of the spin orbit coupling, and hence 
F z . Both analytical and numerical results are in excellent 
agreement with each other. This result survives at an 
arbitrary magnitude of the Berry curvature. 




FIG. 3: Distribution of particles over the momentum space 
in the stationary state with electric field applied along the 
y— axes at zero Berry curvature. Parameters are as follows. 



kf 



1. E v 



0.003, E x = 0, F z 



0. Average distance 
h = e = 1. 



between impurities is I ~ 0.5 in units m 



In FigfJ] and Fig0] we visualized the stationary dis- 
tribution function of particles in the momentum space, 
calculated numerically for the system placed in the uni- 
form electric field. Brighter areas represent higher den- 
sities. To increase the contrast, we subtracted the initial 
equilibrium distribution (i.e. without electric field) from 
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equations (JTJ, (J2J. The new system lf2*7|) already 

corresponds to the classical Hamiltonian evolution with 
the following Hamiltonian 



-1 



o 



FIG. 4: Distribution of particles over the momentum space 
in the stationary state with electric held applied along the 
y— axes at the non-zero Berry curvature. Parameters are as 



follows, ki 



1, E v 



0.003, E x 



0, F z 



distance between impurities is I ~ 0.5. 



-1. Average 



the calculated one. In Fig|2| the Berry curvature is set 
to zero. In this case the correction to the equilibrium 
distribution is symmetric along the electric field. Fig0| 
shows that when the Berry curvature becomes nonzero, 
the distribution function acquires the additional asym- 
metric contribution, which obviously deposits to the Hall 
current. 



E. Analogy between extrinsic and intrinsic AHE 

As we showed above, in the AHE induced by Berry 
curvature there is a possibility of disorder effects similar 
to previously studied extrinsic contributions to Hall cur- 
rent due to gradient terms of the disorder potential. In 
this subsection we discuss how this relation is encoded 
in the wave packet equations QJ, (J5J). We consider here 
only the case of a constant Berry curvature. The wave 
packet dynamic equations do not arise from a particular 
Hamiltonian in a sense that coordinates and momentums 
of a wave packet are not canonical variables. Neverthe- 
less, it is possible to make them arise from a Hamiltonian 
by adding terms of higher order in the potential gradient. 
We remind the reader that equations and 10) were de- 
rived in the approximation of a smooth potential so that 
its higher gradients could be discarded. Therefore such a 
procedure should not change the physics at least to the 
order of the approximation of the whole theory. Consider 
the following system 

-^k = eE - W(r) - V(k x W(r))F (26) 
at 



H(k, r) = e(k) - F(k x (eE - W(r)) + V(r) - eEr (28) 

This Hamiltonian has the same structure as in typical 
models of the extrinsic AHEi^i^i^ and the Berry curva- 
ture now plays the role of an effective spin orbit coupling 
constant. As we mentioned, the main contributions to 
the extrinsic Hall current have been proved to be the 
side jump and the skew scattering. In our previous cal- 
culations, including numerics, we ignored the skew scat- 
tering mechanism, nevertheless the skew scattering can 
be calculated by same weave packet techniques because 
the second term in the rhs. of l|28|) is responsible for 
both skew scattering and the side jump contributions in 
the extrinsic AHE. Note also that the skew scattering 
current will depend on the Berry curvature and thus will 
also have a geometric interpretation. The fact that the 
side jump is the only other surviving contribution in the 
extrinsic AHE is in agreement with our previous findings. 



IV. RASHBA 2DEG WITH EXCHANGE 
COUPLING 

Recently studied systems with AHE such as DMS and 
Rashba 2DEG have Berry curvature that strongly de- 
pends on momentum of the wave packed. When F z 
is no longer a constant the simple arguments leading to 
canceling of some terms may not work. One can notice 
that when an electron decelerates it has in general differ- 
ent momentum than when it accelerates, in other words 
the uniform electric field accelerates all electrons down 
to the bottom of the Fermi sea but impurities produce 
nonzero deceleration in average only for electrons near 
the Fermi surface. The Berry phase acquired by acceler- 
ating electrons depends now not only on the acceleration 
but on the momentum itself therefore contributions from 
impurities and from uniform electric field not necessarily 
cancel each other. 

The Hamiltonian of R2DEG with the electric field 
along the x-direction and the exchange field in the z- 
direction is 



H = H a + H, 



so 



H, 



exch 



- eE x x + Vi 



imp 



bk 2 

H = — cr , b = 1/m 



(29) 



(30) 



<9e(k) 



(eE - W(r)) x F 



(27) 



dt dk 

Except for one term of the second order in gradients of 
the potential the system is equivalent to the wave packet 



Hso = ^(k y a x - k x a y ) 



(31) 



(32) 
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where Vi mp is the impurity potential. 

Diagonalizing the unperturbed part of the Hamiltonian 
we find the energy dispersion 



e±(k) = bk 2 ± \/{Xk) 2 + h 2 



(33) 



There are two bands, the minor band has the plus sign in 
the above expression. Fermi momentums can be derived 
by inversion of the equations 133|) . We denote them as 
fc F _ and fc F+ respectively for the major and minor bands. 
The Berry curvature is different for different bands 



F! = 21m 



du s 
dk v 



du s 
dk x 



(34) 



here \u s ) is the Bloch state in the absence of electric field 
and impurities and s = ± is the index of the band: the 
plus is for the minor and the minus is for the major one. 
In the case of the Hamiltonian (|29|) the Berry curvature 



t ,30 



(see also Appendix for an alternative derivation) 



hX 2 



2[(Xk) 2 + h 2 } 3 / 2 



(35) 



at Afc << h this gives F^(k) rs —sX 2 /(2h 2 ) — const and 
we reduce the problem to the one considered in previous 
sections. In the opposite limit at Afc >> h we find 



2Xk 3 



(36) 



The anomalous velocity has the same magnitude but 
different signs for different bands so filled states with the 
same k from different bands will compensate each other 
in calculations of the clean contribution and to find the 
Hall current we should only integrate over uncompen- 
sated states of the major band. The Hall current from the 
unperturbed electron distribution function in the clean 
case is 



27T k f 



■(clean) _ 
Jyx 



kdk 
J2^T 2 



(eE x F-(k)) 



k 



2-rr k F ~ 



kdk eE x hX 2 
(2ir) 2 2[(Afc) 2 + h 2 } 3 ' 2 

< J E i I I 



k F + 
,2 l 



2(2tt) \ ^/l + (k F+ X/h) 2 



1 



^l + {k F -X/h) 2 



(37) 



in the case Ep » Afc F ± >> h the expression l|37|) is 
simplified 



■(clean) 
Jyx 



\/ ^ ip k F \ e 2 E x h(k F ^ - k FA 

AK <C Hip "~ — 



2m 



(27r)2Afc 2 



{(fc F _-fc F+ )«2Am} = ^^ (38) 



note that here kp + is the Fermi momentum of the mi- 
nor band. The above formulas are valid when there are 
electrons in both bands. At a sufficiently low Fermi level 
(Ep < h,Xk F -) the minor band becomes depleted and 
the transport is only in the major band. In this case the 
clean limit Hall current is 



■ (clean) 
Jyx 



e 2 E x 
2(2vr) 



f 



v/1 + (k F .X/h) 2 



(39) 



We address next the effects of impurity scattering within 
the formalism that we have developed. 

As in the case of the constant Berry curvature, the 
non-equilibrium contribution to the distribution function 
is the sum of the usual one g s (k) , which appears to com- 
pensate the electron acceleration 



.g s (k)=e^r||(-^)«i~° 



(40) 



and of the anomalous one g s p, (k) due to the shift of 
the kinetic energy of the scattered particle. We sup- 
pose here that the impurity potential is weak, namely 
V(r) << E F ,Xk F and disregard the variation of the 
Berry curvature near a wave packet during a scattering 
because the absolute magnitude of the momentum does 
not change appreciably. The derivation of the anomalous 
correction is analogous to (|I0|I - 1I3|) . 



g%(k) = -JkeE x k y Fl 
The normal velocity of the wave packet is 



v 



(normal) 



de s (k) 
dk x 



(41) 



(42) 



The expression for the drift velocity due to the side 
jumps along the y-axes is the same as in (|16|) . The side 
jump current is 



M) - 

Jyx 



£ e |<„«>(k) W < k )*g^ 



(k)k 2 S{E f 



e.fc 



he 2 E r X 2 



k 2 



2(4tt) \[(Xk F+ ) 2 + h 2 } 3 / 2 



kp_ 



[(Afc F _) 2 + /i 2 ] 3 / 2 



(43) 



At Ep >> Xkp >> h this expression can be simplified 

he 2 E T 



Jyx 



8irE F 

If Ep < h only the major band survives and 
-he 2 E r X 2 



Jyx 



kp_ 



2(4tt) [(Afc F _) 2 + h 2 } 3 / 2 



(44) 



(45) 
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FIG. 5: (Color online) Hall conductivity vs the Fermi energy 
in the R2DEG. Black line is the pure limit Berry phase pre- 
diction. The red curve shows the prediction of the theory 
that includes impurity scatterings (1431 . 14711 . Data are given 
in units e = h = m = 1 for the case h = 0.7, A = 0.4. 
Sharp change in the behavior above Ef = 0.7 is due to the 
appearance of electrons in the minor band. 



The normal current contribution is 



■ (normal) 
Jyx 



d*k (normal) s 



5|(k) 



= E s=± / / de sM e 2 E x F!{k)kH{E F - e,, k ) 

(46) 

Comparing the last expression with the one for the side 
jump current (|43|l we find that they are different only by 
an exchange of cos(<yj) and sin(<^) under the integral over 
the angle. Thus we arrive at the general result 



•(normal) 
Jyx 



Jyx 



Finally, 



■(total) _ -(clean) , -(sj) , • 
Jyx Jyx ~ Jyx 'J 



(norma 
yx 







(47) 



(48) 



As it is seen from 11440 . 1)380 the side jump contribution 
may not change the sign in comparison with the clean 
contribution for sufficiently large Fermi energy. This is 
very distinct from the case of the constant Berry curva- 
ture where such a contribution exactly cancels the pure 
one. At the limit Ep >> Xkp >> h we find for the total 
current 



■(total) 
Jyx 



he 2 E x 
2ttE f 



(49) 



The total Hall current not only has the same sign, as the 
clean limit prediction, but also increases due to scatter- 
ings. Note also, that up to a coefficient of the order unity, 
in this limit, the total current has the same dependence 
on various parameters. 

In Fig|3]we compare the theoretical prediction of Eqs. 
(|48|l with the prediction of the pure Berry phase contri- 
bution {57)1 . 



V. DISCUSSION AND CONCLUSION 

We calculated the anomalous Hall current in the pres- 
ence of a smooth disorder potential by combining the 
wave packet dynamic equations with the classical Boltz- 
mann equation. For the case of a constant Berry curva- 
ture our results are in agreement with some known pre- 
dictions of the very earliest works on AHE. The gauge 
invariance of our approach clarifies the physical meaning 
of various disorder contributions to the Hall current and 
reveales the relation between the extrinsic, related to the 
gradients of the disorder potential, and Berry curvature 
effects in the presence of disorder. 

We have shown that early results, like that of the Lut- 
tinger's sign reversal prediction, cannot be used directly 
for new applications such as diluted magnetic semicon- 
ductors, where the Berry curvature is strongly momen- 
tum dependent. 

Our prediction for the Hall current FigEl is distinct 
from the clean limit (or the pure Berry phase) one. Thus 
the total Hall current changes the sign when the Fermi 
level increases. In contrast, the clean limit prediction 
for the models considered remains of the same sign. The 
filled states in both bands with the same momentum com- 
pensate each other, hence for the clean case the band 
having more filled states always wins. In contrast, impu- 
rity contributions are sensitive not only to the density of 
states near the Fermi level, but also to the magnitude of 
the Berry curvature near the Fermi level. 

The change of sign of the Hall current as a function 
of chemical potential originates from the competition be- 
tween two bands. The minor band has the smaller num- 
ber of filled states than the major one but it also has 
lower Uf and hence the stronger F z (kp). According to 
our calculations, such a change of sign should generally 
happen when the chemical potential is close to the deple- 
tion point of the minor band. When the minor band is 
totally depleted we always observe the reversal of the sign 
of the Hall current in comparison to the clean limit con- 
tribution. We also note that such a competition should 
be applicable to the skew scattering mechanism as well. 
This follows from the fact that the strength of the skew 
scattering should be proportional to the Berry curvature 
at the Fermi surface. 

Our results can be generalized to DMS and although 
in this case the complicated behaviour of the bands, as 
compared to the R2DEG, may change the quantitative 
behavior there is still a competition in DMS of Hall cur- 
rents from major and minor bands and Berry curvature 
is also usually decreasing when momentum is increasing. 
However the intricate dependence of the clean limit AHE 
observed in these systems versus the simple dependence 
observed in the R2DEG will require a careful analysis. 

Although we have verified most our predictions by 
numerical simulations, the correspondence between the 
wave packet approach and the evolution of a true quan- 
tum system with disorder remains to be investigated. An 
effect not taken in to account in our calculations is that 
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the spin density in the wave packet polarizes when the 
wave packet is accelerated^. This effect is crucial for 
the Spin Hall Effect and its importance for calculations 
by the Boltzmann equation should be understood. 

The final theory of disorder in the AHE can be achieved 
only by purely quantum mechanical calculations for ex- 
ample, based on the Keldysh technique. At present 
new approaches to systems with Berry curvature have 
emerged^, as well as new mechanisms to generate the 
anomalous Hall current^, but we believe that the wave 
packet approach combined with the classical Boltzmann 
equation is worth studying because it provides the sim- 
plest, to our knowledge, demonstration of the related 
physics and the wave packet equations themselves have 
been well justified by the quantum theory 
Note added. After completion of this work we became 
aware of the related effort by V. K. Dugaev et al£^ of a 
more quantum mechanical approach to the problem with 
similar qualitative conclusions. Our approach however is 
quite different and in combination with this work may 
shed further light on the physical interpretation of the 
effect. 
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APPENDIX A: CLASSICAL INTERPRETATION 
OF THE WAVE PACKET DYNAMICS IN A 
RASHBA 2DEG 

To make the semiclassical result more transparent we 
show how the anomalous velocity arising from the Berry 
curvature can appear in a purely classical system. This 
classical anomalous velocity originates from the non- 
adiabatic contributions to the equations of motion in the 
linear response regime and are not present when a sim- 
pler adiabatic approximation is considered. 

We consider the motion of a classical particle having 
the electric charge e and the classical spin S attached to 
it. Choose the Hamiltonian to be 

k 2 

H = h 2X(k y S x - k x S y ) - eE x x + 2hS z (Al) 

The analogy with the Rashba Hamiltonian (|29|l is obvi- 
ous. In the classical problem we substitute operators by 
corresponding classical variables. We will suppose also 
that \S\ = 1/2 as for true electrons. The Hamilton equa- 
tions for the evolution of coordinates (x, y) and momen- 
tum (k Xl k v ) and Landau-Lifshitz equations for the evo- 
lution of spin components read 

x=—- 2XS V (A2) 



y = ^- + 2XS X (A3) 

TO 

k x (t)=eE x (A4) 

k v {t) = (A5) 

S = -A x S (A6) 



where A is the effective magnetic field acting on the spin. 
From the Hamiltonian (|A1|) we can read this effective 
magnetic field acting on the particle with the momentum 
k = (k x , ky) 

A = —2{Xk y , ~\k x , h) (A7) 

The absolute magnitude of this field is 

|A| = v/4(/c 2 A 2 + h?) (A8) 

Assuming that the electric field E x is weak, the varia- 
tion of the momentum and hence of the effective magnetic 
field is very slow and, if at the initial moment the classical 
spin is directed along the magnetic field, it will follow the 
direction of this slowly varying field. This adiabatic ap- 
proximation follows from the Landau-Lifshitz equations 
and yields an approximate solution for (|A6|) given by 

S^(t) = ±|S|£®, where "+" and "-" correspond to 
the initial direction along or opposite to the magnetic 
field respectively. Substituting this time dependence of 
the spin direction into equations for particle velocities 
one can find that in the adiabatic approximation 



(A9) 



where eo(k) = j^Ty/{k\) 2 + h 2 , i.e. a particle with spin 
along the field A has energy dispersion as an electron 
in the major band of R2DEG. In the strict adiabatic 
approximation of this classical Hamiltonian particles do 
not exhibit the anomalous velocity seen in the second 
term of the group velocity of the quantum-mechanical 
wavepacket dynamics described by Eq. (2). 

However, this classical model does capture the correct 
contribution to the anomalous Hall effect related to the 
Berry curvature once we go beyond the strict adiabatic 
approximation and consider instead the more general lin- 
ear response. The nonadiabatic corrections describe de- 
flection of the classical spin from the direction of the 
instant magnetic field. They are small but can still be 
of the first order in the electric field and hence affect the 
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transport properties in the linear response. Looking for 
a solution of the Landau-Lifshitz equations (|A6I) in the 
form S = S^ a ' + 8S we find that in the first order in 
the external electric field the spin direction acquires the 
component perpendicular both to the effective field and 
its derivative^ i.e. 



SS y ((k) = T \S\ 



\hk„. 



2[(fcA) 2 + /i 2 ] 3 / 2 



(A12) 



If h =/= this correction has in-plane components and 
hence affects the velocity of the particle, i.e. it gives a 
contribution to the right hand side of Eq. IA2I and IA3I 
Substituting ||a7| . l|Al| into (|A10|1 we find 

as -« k > " ^[(W + wv. < A11 » 



Then substituting ijATTJl and (|A"l"2)) into (|X2|) . llAHl for 

\S\ = 1/2 we find equations identical to the wave packet 
ones with the Berry curvature given in Eqs. J5J and l|35(l 

We emphasize that although being an interesting 
demonstration of the effect, the model of classical charged 
particles with classical spins cannot always be employed 
in calculations of the Berry curvature and quantum cor- 
rections for other types of the Hamiltonian may appear. 34 
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